arXiv: 1509.05271 vl [astro-ph.GA] 17 Sep 2015 


Astronomy Astrophysics manuscript no. eeuielion ©ESO 2015 

September 18, 2015 


The Gaia-ESO Survey: 

New constraints on the Galactic disc velocity dispersion and 

its chemical dependencies 

G. Guiglion^, A. Recio-Blanco^, P. de Laverny^, G. Kordopatis^, V. S. Mikolaitis^, I. Minchev^, G. 

Ghiappini^, R. F. G. Wyse^, G. Gilmore^, S. Randich^, S. Feltzing®, T. Bensby®, E. Flaccomio^, S. E. 
Koposov'^’® E. Pancino®, A. Bayo^*^, M. T. Gostado^^, E. Franciosini®, A. Hourihane'^, P. Jofre^, G. Lardo^^, 
J. Lewis^, K. Lind^®, L. Magrini®, L. Morbidelli®, G. G. Sacco®, G. Ruchti®, G. G. Worley^, S. Zaggia^"^ 

^ Laboratoire Lagrange, Universite Cote d’Azur, Observatoire de la Cote d’Azur, CNRS, Blvd de I’Observatoire, CS 
34229, 06304 Nice cedex 4, France 

^ Leibniz-Institut fiir Astrophysik Potsdam (AIP) An der Sternwarte 16, D-14482 Potsdam, Germany 
^ Physics and Astronomy Department, Johns Hopkins University, 3400 North Charles Street, Baltimore, MD 21218, 
USA 

Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 OHA, United Kingdom 
® INAF - Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125, Florence, Italy 
® Lund Observatory, Department of Astronomy and Theoretical Physics, Box 43, SE-221 00 Lund, Sweden 
^ INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134, Palermo, Italy 
® Moscow MV Lomonosov State University, Sternberg Astronomical Institute, Moscow 119992, Russia 
® INAF - Osservatorio Astronomico di Bologna, via Ranzani 1, 40127, Bologna, Italy ASI Science Data Center, Via 
del Politecnico SNC, 00133 Roma, Italy 

Instituto de Flsica y Astronomha, Universidad de Valparailso, Chile 
Instituto de Astroflsica de Andalucia-CSIC, Apdo. 3004, 18080, Granada, Spain 

Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, United 
Kingdom 

Department of Physics and Astronomy, Uppsala University, Box 516, SE-751 20 Uppsala, Sweden 
INAF - Padova Observatory, Vicolo dell’Osservatorio 5, 35122 Padova, Italy 

Received 13/02/2015; accepted 26/08/2015 


ABSTRACT 

Context. Understanding the history and the evolution of the Milky Way is one of the main goals of modern astrophysics. 
In particular, the formation of the Galactic disc is a key problem of Galactic archaeology. 

Aims. We study the velocity dispersion behaviour of Galactic disc stars as a function of the [Mg/Fe] ratio, which for 
small metallicity bins can be used as a proxy of relative age. This key relation is essential to constrain the formation 
mechanisms of the disc stellar populations as well as the cooling and settling processes. 

Methods. We used the recommended parameters and chemical abundances of 7 800 FGK Milky Way field stars from 
the second internal data release of the Gaia-ESO spectroscopic Survey. These stars were observed with the GIRAFFE 
spectrograph (HRIO and HR21 setups), and cover a large spatial volume in the intervals 6 < 7? < 10 kpc and \Z\ < 2 kpc. 
Based on a chemical criterion, we separated the thin- from the thick-disc sequence in the [Mg/Fe] v.s [Fe/H] plane. 
Results. From analysing the Galactocentric velocity of the stars for the thin disc, we find a weak positive correlation 
between Vtf, and [Fe/H] that is due to a slowly rotating [Fe/H]-poor tail. For the thick disc stars, a strong correlation 
with [Fe/H] and [Mg/Fe] is established. In addition, we have detected an inversion of the velocity dispersion trends 
with [Mg/Fe] for thick-disc stars with [Fe/H] < —O.lOdex and [Mg/Fe] > -1-0.20 dex for the radial component. First, 
the velocity dispersion increases with [Mg/Fe] at all [Fe/H] ratios for the thin-disc stars, and then it decreases for the 
thick-disc population at the highest [Mg/Fe] abundances. Similar trends are observed for several bins of [Mg/Fe] within 
the errors for the azimuthal velocity dispersion, while a continuous increase with [Mg/Fe] is observed for the vertical 
velocity dispersion. The velocity dispersion decrease agrees with previous measurements of the RAVE survey, although 
it is observed here for a greater metallicity interval and a larger spatial volume. 

Conclusions. Thanks to the Gaia-ESO Survey data, we confirm the existence of [Mg/Fe]-rich thick-disc stars with cool 
kinematics in the generally turbulent context of the primitive Galactic disc. This is discussed in the framework of the 
different disc formation and evolution scenarios. 

Key words. Galaxy: abundances - Galaxy: disc - Galaxy: kinematics and dynamics - Galaxy: stellar content - Stars: 
abundances 
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1. Introduction 

During the past decades, strong observational efforts have 
been made to understand the formation and evolution of the 
Milky Way (MW). To identify and characterise the Galactic 
structures, we need to study two main physical properties: 
the stellar chemistry and the kinematics. The photospheric 
elemental composition of low-mass stars reflects the chem¬ 
ical evolution history of the stellar populations, while the 
3D velocities are the fossil records of the stellar motions 
inside the Galaxy. Thus, massive spectroscopic surveys are 
the best tools for distinguishing stellar populations. 

The first massive and low-resolution surveys such as 


RAVE (Steinmetz et al. 


and SEGUE (Yanny et al. 


20091 provided pioneer results with excellent statistics, but 


more detailed chemical information can be obtained from 
survey s with a higher reso lution, such as the Gaia-ESO 
(GES , Gilmore et al. 20121 or APOGEE (Allende Prieto 
et al.|^08p. 


posed by several studies ( 

Hernquist & Quinn [1989 

Quinn 

et al. 1993 

Walker et al. 

1996 

Villalobos & Helmi 

2008p. 

Abadi et al 

(2003P argued that 

the thick disc is composed 


while the coalescence of a gas-ric h satellite forming a thic k 
disc in situ was also suggested by [Brook et al. (2004 2007 k 


More recently, radial migration ([Sellwood & Binney 
2002| , driven by the co-rotation (churning) or the Lindblad 
(blurring) resonances with the disc spiral structure, has 
been evoked to explain the eme rgence of the thick disc. Nev¬ 
ertheless, several authors (e.g. Minchev et al. ][M^[MartIg 
et al. 2014 , jVera-Giro et al. 2014p used numerical models 
to propose that migration in the sense of churning does not 
significantly contribute to disc thickening. 


the Geneva-Gopenhagen Survey, Nordstrom et al. (20041 
found that the AVD relation of the young disc (younger 
than 7 Gyr) follows a power law. This was later revised by 


Seabroke & Gilmore (2007), however, who showed that no 
power law is needed to explain the shape of the relation. Fi¬ 


nally, Sharma et al.|(2014l constrained the thin-disc AVD 


relation as a power law with RAVE data by fitting a kine¬ 
matic model to the observations. 

All these studies are based on small samples with well- 
constrained distances (Hipparcos volume) or on larger sam¬ 
ples with a low or intermediate spectral resolution. How¬ 
ever, determining precise ages for a larger star sample is 
quite challenging because we lack good distance measure¬ 
ments and precise fundamental parameter determinations. 
For the youngest and oldest stars, moreover, the isochrones 
are crowded in these regimes, and it is harder to derive pre- 


The existence of two se parated thin and thick discs is 
still a matter of debate (see Bovy et al. |2012a| ) in a context 


gested as an age proxy by several authors ( 

Matteucci|2001 

Bovy et al. 

2012c Minchev et al. |2014b) 

Indeed, in the 


where several scenarios are proposed to explain the forma¬ 
tion of the Galactic disc stellar populations. The vertical 


ejecta rapidly enriched the interstellar medium and fixed 
the [a/Fe] ratio at a plateau value. Because they evolve on 
a long timescale, SNIa generate more iron, which leads to 
a decline of the [a/Fe] ratio with increasing [Fe/H] and 


forms a typ ical b reak (or knee). In t eresti ngly, Haywood 
et al. (2013) and Bergemann et al. (2014|) found a tight 


On the other hand, Bournaud et al. (2009) showed that 
stellar scattering by massive clumps could explain the cre¬ 
ation of the t hick disc through a pur ely internal mecha¬ 
nism. Finally, Haywood et al. (2013) proposed that the 
early epoch of the disc was driv en by an inten s e star for- 
mation, which was q uantified by Lehnert et al. (2014) and 
Snaith et al. (2015). All these proposed scenarios can be 
constrained by studying the chemo-dynamical relations fol¬ 
lowed by the stellar populations of the disc in a large spatial 
volume and with robust statistics. 

To understand the dynamical heating of disc stars, one 
can investigate the age-Velocity-Dispersion (AVD) relation 
in the very close Solar neigh bourhood, which was first stud¬ 
ied by Wielen ( 1974| 1975), who used one thousand stars 
within 20 pc from tne Sun. These works highlighted the 
fact that o ld stars have a higher velocity dispersion. jWie-j 
len (1977) showed that this interesting behaviour could 
be explained by the diffusion of the orbits due to interac- 
tions with the l o cal flu ctuations of the gravitational field. 
Garlberg et al. (1985) provided a new AVD relation for 
about 500 h'-type stars, showing a fiat behaviour with age 
for stars older than 6 Gyr. Based on the dwarf sample of 


correlation b etween the stellar ages and the [Mg/Fe] ratio. 
We note that Lee et al. (2011a) observed a steeper increase 
of the velocity dispersion as a function of the [Mg/Fe] ratio 
in the thick disc than in the thin component. 

More recently, a relation between the velocity disper¬ 
sion and the [Mg/Fe] ratio, used as an age proxy, has been 
presented by Minchev et al. ( 2014b[) . They derived a chemo- 
kinematical relation for the Galactic disc for about5 000 gi¬ 
ant stars of the Solar suburb (|^| < 0.6 kpc and 7 < R < 
9kpc) from the RAVE DR4 survey ( Kordopatis et al.|2013 ). 
They found a surprising behaviour for the most Fe-poor and 
Mg-rich stars ([Fe/H] < —O.Sdex and [Mg/Fe] > -|-0.4dex): 
a strong decrease of the velocity dispersion is observed while 
it is expected to increase because these stars should corre¬ 
spond to the oldest population of the Galactic disc. We 
note that they also observed this decrease in the SE GUE 
G-dwarf sample (R = 2000) from Lee et al. (2011a) and 
the [a/Fe] ratios, instead of [Mg/Fe]. In our study, we fo¬ 
cus on the results obtained with RAVE because they have a 
higher resolution and were obtained using the [Mg/Fe] ra- 
tio. By comparing their results with numerical simulations. 


* Based on observations collected with the FLAMES spectro¬ 
graph at the VLT/UT2 telescope (Paranal Observatory, ESO, 
Chile), for the Gaia-ESO Large Public Survey, programme 
188.B-3002 


Minchev et al. (2014b) concluded that this unexpected be- 
haviour is due to perturbations by massive mergers on the 
outer part of the disc in the early epoch of our Galaxy, trig¬ 
gering a strong radial migration process that redistributed 
stars form the inner disc into the Solar suburb on near¬ 
circular orbits and lower velocity dispersion. 

To extend this last study, we propose here to analyse the 
behaviour of the velocity dispersion of the Galactic disc as 
a function of its [Mg/Fe] ratio with GES data. The observa¬ 
tions were made with the VLT/UT2 of the European South¬ 
ern Observatory. Our analysis is based on the data from the 
second internal data release (iDR2), which is composed of 
about 8 000 FGK dwarf and giant stars. We here have the 
opportunity of probing the Galactic disc more deeply be¬ 
cause GES observes fainter targets farther away from the 
Solar neighbourhood. Moreover, we have a higher precision 
on the atmospheric parameters and chemical abundances 
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thanks to the higher spectral resolution and wavelength 
coverage than was possible with RAVE. 

The paper is organised as follows: we first present our 
stellar sample and chemically define the Galactic discs 
(Sect.|^, then we detail our method for the velocity disper¬ 
sion derivation (Sect.[^. Section|^is devoted to presenting 
our results, and we compare our chemo-kinematical rela¬ 
tion with the one of the RAVE survey in Sect. Finally, 
we discuss the constraints provided by this study on the 
different scenarios of the Galactic disc formation in Sect. |6] 
and then conclude and summarise our work in Sect. 0 

2. Stellar sample and characterisation of the 

Galactic disc 



Our work is based on the second internal data release 


iDR2) of the Gaia-ESO spectroscopic Survey (Gilmore 
et al. 2012|). We chose a sample of FGK Milky Way held 
stars observed by the GIRAFFE spectrograph in both 
HRIO (R ~ 19 800) and HR21 (R ~ 16 200) setups. In par¬ 
ticular, we selected stars with available radial velocity Vadi 
effective temperature T^a, surface gravity 1 ^( 5 ), global 
metallicity [M/hQ and [Fe/H] and [Mg/Fe]j abundance 
ratios; this results in a sample of 7 800 stars. 

We recall that for iDR2, the radial velocity Vad was 
determined in a first step (complete procedure described 
in Koposov et al., in preparation), with a median error of 
about 0.3km s“^. Then, the fundamental parameters Teffi 
log(g), [M/H] and [a/Fe] were computed in three d ifferent 
independent proced ures (see a short description in Recio- 
Blanco et al.[ 2014). The median dispersions among these 


independent measures for our sample are 41K, 0.09 dex, 
and 0.05dex for T^s, log( 5 ), and [M/H], respectively. We 
selected stars with a minimum signal-to-noise ratio (S/N) 
equal to 10, leading to distributions with a median value of 
23 in HRIO and 50 in HR21. 

From the atmospheric parameters, a and iron-peak el¬ 
emental abundances were computed by combining the de¬ 
terminations of three distinct methods (see Recio-Blanco et 
al., in preparation). To chemically characterise the Galactic 
disc, we focus in the following on the recommended mag¬ 
nesium and iron abundances. We adopted the definition 
[Mg/Fe] = [Mgi/Fei], rejecting the Fell abundances that 
are measured from a single spectral line. In the adopted 
sample, the typical median dispersions of the three indepen¬ 
dent measures are 0.06 and 0.04dex for [Fe/H] and [Mg/Fe], 
respectively. 

To determine the stellar kinematics , we f ollow ed the 
procedure de s cribed in Kordopatis et al. (20111 and Recio- 
Blanco et al. (2014|). We recall that we derived the absolute 
magnitudes and then the line-of-sight distances, projecting 
Teff, log(fl) and [M/H ] on a set of Yonsei-Yale isochrones 
( Demarque et al.|2004| ) . The three-dimensional Galactic co- 
ordinates of the stars were then calculated. Finally, the 
Galactocentric radial, azimuthal, and vertical velocities 
(Vfl, V^, and Vz, respectively) were derived for all the 
stars with available proper motions from the PPMXL cata¬ 


logue (Roeser et al. 


of Rest (LbR) (U, 


2010). We adopted the Local Standard 
W)q = (11.1,12.24, 7.25) km s-^ 


|Mg/F.l-log[5@ff]_^- log[:^] 


Fig. 1. [Mg/Fe] as a function of [Fe/H] for the 1 230 iDR2 stars 
with S/N > 40 (16% of the selected sample). The error bars 
show the typical median dispersion. The red solid line represents 
our division between the thin- (below) and the thick- (above) 
disc stars. 


Schbnrich et al. ( 2010| . The LSR is presumed to be on a 
circular orbit with an azimuthal velocity V/ = 220 km s“^. 
For the two components of the proper motion and fxg, 
the typical errors are about 7 mas/yr. The resulting median 
errors on the whole sample are 52, 39 and 52 km s“^ for V/j, 
V^, and V^, respectively. Finally, our main sample is com¬ 
posed of 6 800 stars with fundamental parameters, [Mg/Fe] 
abundances, and 3D kinematics. 

We have explored the existence of a chemical gap be¬ 
tween the thin and the thick disc populations in the 
[Mg/Fe] us. [Fe/H] plane. This gap was first detected in 
high spectral resolution samples in the close Solar neigh¬ 


bourhood by several authors (e.g. Fuhrmann 1998[ [Feltz- 
~ ' ' has been confirmed 


ing et al(]|2Cioi| Bensby et al. 112003 ). it has b( 
more recently by Adibekyan et al. ( 2011| and 


(2014) still in the close Solar n'erghbourhood and beyond 


Bensby et al. 


thanks to GES iDRl (Recio-Blanco et al. 2014 Mikolaitis 
et al. 2014) and APOGEE (Nidever et al. 2014) data. To 
define this gap, we follo wed the same procedur e as in|Recion 


of 


Blanco et al. (2014) and Mikolaitis et al. (2014). We selected 
every star with S/N > 40, (1 230 stars corresponding to 
16% of the selected sample) that we decomposed into nine 
[Fe/H] bins from -1.50 to -1-0.50 dex. Then, we searched the 
minimum count of each bin in the corresponding [Mg/Fe] 
distribution. We obtained a quasi-linear separation from 
[Fe/H] = —1.25 to —0.1 dex. For the [Fe/H]-rich and -poor 
end, no gap has been detected, and we thus assumed a 
plateau (see Fig.[^. In the following, we identify stars with 
[Mg/Fe] below that separation as thin disc stars and those 
above that separation as thick disc stars; we rejected stars 
with [Fe/H] < —1.25 dex. The possible data contamination 
by halo stars is discussed at the end of Sec t. |4.3| The sep¬ 
aration foun d by Recio-Blanco et al. (2014) and |Mikolaitis 
(201^ with GES iDRl data is thus again confirmed 
E; 


et al. 
with the G 


^ [M/H] corresponds to the proportion of all elements heavier 
than helium 


O 


S iDR2 data. We emphasize the fact that the 
two iDRl thin- to thick-disc separations were found in the 
[a/Fe] — [M/H] and [Mg/M] — [M/H] plane, respectively, 
where M is the global metallicity. As a consequence, our 
[Fe/H] — [Mg/Fe] division is slightly different from these 
two previous studies, but totally compatible. 
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3. Deriving the velocity dispersion of the Galactic 

disc 

This section describes in detail our estimation of the stellar 
velocity dispersion as a function of the chemical abundance 
[Mg/Fe], which was adopted as a stellar clock. 

For the main sample, we decomposed the [Fe/H] — 
[Mg/Fe] plane into 2D bins. First, the [Fe/H] axis was cut 
from [Fe/H] = —0.65 to [Fe/H] = +0.15 dex into four bins 
with a constant 0.2 dex step. The [Fe/H]-poor less popu¬ 
lated tail was included in a single bin from [Fe/H] = —1.25 
to —0.65 dex, and we adopted the same procedure for the 
[Fe/H]-rich end from [Fe/H] = +0.15 to +0.65 dex. In to¬ 
tal, six [Fe/H] bins were thus considered. Second, we cut the 
[Mg/Fe] axis into bins of 0.1 dex starting from [Mg/Fe] = 
—0.2 dex up to [Mg/Fe] = +0.4 dex, and we considered a 
last bin including all the stars with [Mg/Fe] > +0.4 dex. 
We note that for each [Fe/H] bin, the thin- to thick-disc 
transition presented in Sect, [^divides one [Mg/Fe] bin into 
two distinct sub-samples. 

For a given 2D bin, the average Galactocentric veloc¬ 
ity fi and the intrinsic velocity dispersion tJi were derived 
by model fitting with a maximum likelihood approach. We 
computed the likelihood as follows: 


N 

2=1 


1 

v/27r(crf + evf) 


exp 


2(af + evf)J’ 


( 1 ) 


where Vi is the measured Galactocentric velocity for a 
star in a given velocity component and evi its measured 
error. The measured velocity distribution can be seen as 
the convolution between the intrinsic distribution and the 
measured errors. In an equivalent manner, we minimised 
the log-likelihood A = — 21nL with respect to both the 
mean velocity p and the intrinsic velocity dispersion ai. 
To find the minimum, we investigated where the first-order 
derivative of A goes to zero, which simultaneously solves 
the equations f^= 0 and 0 , resulting in 


N 

E 


N 


-a^E 




= 0 , 


( 2 ) 


and 


N 

E 


(a} + - {vj - pY 


= 0 , 


( 3 ) 


respectively (see also for instance Godwin & Lynden- 


Bell 1987 Pryor & Meylan 19931. This was performed using 
IDL procedures. 

The quality of the kinematical data is crucial for the 
velocity dispersion estimation because the errors in the dis¬ 
tances and proper motions propagate in the velocity mea¬ 
surements and sometimes create strong outliers. Therefore, 
we additionally cleaned our main sample so that it includes 
only stars with an accuracy in V/j, , and Vz better than 
30 km s“^. We rejected 5% of the stars that had extreme 
values of the velocity distribution, which were too far from 
the mean velocity. Finally, we estimated the error on ai 
, as the standard deviation of 1 000 bootstrap realisations. 
For a given bin B containing N stars, we created a new bin 
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Fig. 2. Top left: [Mg/Fe] as a function of [Fe/H] for the main 
sample. An example [Fe/H]-rich bin is shown in green (232 
stars). Top-right: Galactocentric radial velocity distribution for 
this bin with Umeas = 49 km s“^. Bottom left: Distribution of 
errors on Vr. Bottom right: Velocity dispersion distribution 
corresponding to 1 000 bootstrap realisations, characterised by 
(Te = 2 km s“^. The intrinsic velocity dispersion derived by max¬ 
imum likelihood is or = 35 km s“^. 


composed of 1 000 x B stars. Then, we randomly selected 
1 000 times N stars in this new bin to measure 1 000 values 
of ai. 

As an illustration, we show in Fig. the steps of the 
procedure we used to retrieve the velocity dispersion in the 
radial component for a given bin in the [Mg/Fe] vs. [Fe/H] 
plane composed of 232 stars with ([Fe/H]) = +0.01 dex and 
([Mg/Fe]) = —0.05 dex. Its velocity distribution is shown in 
the top right panel. We find a measured velocity dispersion 
of (Tmeas = 49 km s“^. By model fitting with the maxi¬ 
mum likelihood approach, we obtain ap = 35km s“^. A 
distribution of 1000 bootstrap realisations is shown in the 
bottom right panel, characterised by a standard deviation 
of (Te = 2 km s“^. 

We have tested the possible influence of the low statis¬ 
tics on the measured velocity dispersion. To this purpose, 
we selected a well-populated [Mg/Fe] interval from the bin 
([Fe/H]) = —0.54dex (azimuthal component) composed of 
100 stars (([Mg/Fe]) = +0.17dex, a^ = 2km s“^). 
We randomly selected 1000 times a percentage of the stars, 
from 5 to 100% with a 5% step. Then, we performed 1000 
derivations of the velocity dispersion, computing an aver¬ 
age and standard deviation. The resulting average of a^ is 
shown in Fig. as a function of the percentage of the se¬ 
lected stars. We see that (cr^) converges to a plateau after 
10% around the reference value a^ = 28 km s“^, with an 
increase of 30 km s“^ for the lowest percentage (composed 
of only five stars). We also tested some other bins at dif¬ 
ferent [Fe/H] and [Mg/Fe] ratios and observed the same 
behaviour. Based on this test, we conclude that bins with 
more than 15 stars are reasonable for a scientific applica- 
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Fig. 3. Average azimuthal velocity dispersion as a function of 
the selected percentage of a bin composed of 100 stars, char¬ 
acterised by ([Fe/H]) = —0.54dex, ([Mg/Fe]) = -|-0.17dex, and 
= 28 ±2 km s“^. The error bars correspond to the standard 
deviation of 1000 bootstrap realisations. 


tion and are suitable for a reliable interpretation in terms 
of kinematical properties of the disc. 

Finally, from our selection criteria in terms of Galac- 
tocentric velocities, we obtained three samples composed 
of 1 998, 2 459, and 2 045 stars for the radial, azimuthal, 
and vertical Galactic velocity component, respectively, as 
shown in Fig. These three samples quite homogeneously 
cover a wide Galactic radius range from 6 to lOkpc, and 
2kpc above and below the Galactic plane. The fraction of 
stars belonging to the thin and thick discs as defined by 
the two sequences in elemental abundances is summarized 
in Table We see that the thin disc represents ~ 80% 
of the sample. We also verified that the mapping of the 
[Mg/Fe] — [Fe/H] plane is the same for the three samples. 

Table 1. Total number of stars (Ntot) for the three velocity 
components {R, (j> and Z), with the respective proportion of the 
thin- and thick-disc samples. 


Velocity component 



Vz 

Ntot 

1998 

2 459 

2 045 

Thin disc 

1593 

1982 

1513 

Thick disc 

405 

477 

532 


4. Velocity gradients and velocity dispersions in 
the Galactic disc 

From the derived velocity dispersions, we first searched for 
correlations between the Galactocentric velocities with re¬ 
spect to [Fe/H] and [Mg/Fe] for the thin and the thick discs. 
We then derived the velocity dispersion relations with re¬ 
spect to [Mg/Fe]. 


4.1. Galactocentric velocities vs. [Fe/H] 

In Fig. 1^ we show the behaviour of the Galactocentric ve¬ 
locities as a function of the [Fe/H] ratio. Each [Fe/H] bin 
is also decomposed into several [Mg/Fe] bins, following the 
same mapping of the [Mg/Fe] — [Fe/H] plane as described 
in the previous section. 


For the radial component {left panel), we see that the 
thick disc shows no trend within the errors with [Fe/H], 
while the thin disc is characterised by a very weak negative 
gradient. We note that the scatter is larger for the most 
[Fe/H]-poor bins in which the number of stars is lower. The 
relation between the vertical Galactocentric velocity and 
the [Fe/H] {right panel) is similar in the thin and thick disc. 
Both thin and thick discs follow a unique sequence with a 
very weak slope, and in addition, the thin disc is compatible 
with no slope. Here again, the scatter of the lowest [Fe/H] 
bins is due to the lower statistics. Finally, the behaviour of 
those two components agrees well with the expected trend 
for the Milky Way thin and thick discs. 


In contrast, the azimuthal velocity {middle panel) shows 
two clearly distinct sequences. The first one corresponds 
to the thin disc with a weak positive gradient equal 
to -1-4 ± 3 km s“^ dex“^. The second sequence, belong¬ 
ing to the thick disc, is clearer with a stronger gradient 
VV^/Vpe/H] = -|-49±10km s“^ dex“^. In both cases, the 
slopes and their associated errors were estimated through a 
linear regression procedure, taking into account the statis¬ 
tical errors of each bin. This azimuthal gradient estimated 
for the thick disc agrees with the one found by |Kordopatis 
et al.| (poTTI (VV^/V[Fe/H] = 45 ± 12 km s"" dex-% 
even though their thick disc w as not based on a chemical 
separation), and more recently, Recio-Blanco et al. (20141 
found VV^/Vpe/H] = -|-43 ± 13 km s“^ dex“^ with GES 
iDRl data. However, the thin-disc azimuthal gradien t dif- 
fers notably from previous studies. Lee et al. (2011bI and 


Recio-Blanco 
—22.6±1.6 km s 


et al. 


(20141 found a negative gradient of 
Hex”^ and —17±6km s“^ 


dex“^, respec¬ 
tively. In our stu dy, the thin disc covers —1 .0 < [Fe/H] < 
-1-0.5 de x, while in Recio-Blanco et al. (2014| and Lee et al. 
( 2011b| ), the range was sniciller, —0.9 <C WH] < 4-0.2 and 
—0.7 < [Fe/H] < -1-0.3, respectively. We first point out that 
the trend with [Fe/H] masks a stronger one with [Mg/Fe], 
which is directly caused by the results being Fe-range de¬ 
pendent. Second, we discuss the robustness of the gradient 
with the Fe-range. Our lowest [Fe/H] bin for the thin disc 
contains 59 stars ('^ 90% of main sequence (MS) stars), 
while the highest [Fe/H] bin contains 299 stars (~ 80% MS). 
We checked that those stars have a chemical pattern in sil¬ 
icon and calcium that is compatible with the one presented 
here in [Mg/Fe]. When we restr ain our thin disc to a [ Fe/H] 
domain similar to that used by Recio-Blanco et al.j (|2014 1 
and Lee et ah] ( 2011bl ), i.e. —0.7 < ]Fe/H] < 4-0.2dex, 
we find a negative gradient equal to —5 ± 5 km s“^ dex“^. 
This result is consistent in terms of slope with the anti¬ 
correlations found previously. In a last check, we decided to 
compute a new gradient, for which we removed the stars in 
the thin-disc sequence with [Mg/Fe] > 0.15 dex. These stars 
are in the [Fe/H]-poor thin- disc tail and show a different 
behaviour (V^ values lower than 210 km s“^, grey and blue 
curves) with respect to the other t hin-disc members. We 
note that Recio-Blanco et al. (20141 did not take this pop¬ 
ulation into account to compute the gradient because they 
lacked good statistics. In this case, we find a slope equal to 
—8±4km s“^ dex“^, which is again consistent with the lit¬ 
erature values. In summary, all these facts indicate that our 
azimuthal gradient is reliable. Therefore, we conclude that 
the thin-disc azimuthal gradient reflects the driving effect 
of [Mg/Fe] because of the dependence on the Fe-range. 
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Fig. 4. Vertical distances from the Galactic plane (Z) as a function of the Galactic radius (R) for the radial, azimuthal, and 
vertical velocity components, composed of 1 998, 2 459, and 2 045 stars , respectively. The black boxes correspond to the spatial 
coverage of the Solar suburb RAVE sample from Minchev et al. (2014b I. We over-plotted the respective normalised distributions 
of R and Z. 
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Fig. 5. Average radial, azimuthal, and vertical velocity (< Vr >, < V 4 , >, and < Vz >) as a function of the [Fe/H] ratio. 
The [Mg/Fe] bins are colour-coded, and the average values of each one are written in the legend. The thin and thick discs are 
represented by open squares (□) and filled circles (•), respectively. The error bars correspond to the standard errors on the mean 
velocity for each bin. The derived thin- and thick-disc gradients are shown as full and dashed lines, respectively. 


4.2. Galactocentric velocities vs. [Mg/ Fe] 

Figure [^presents the Galactocentric velocities as a function 
of [Mg/Fe]. For the radial component {left panel), we found 
no slope for the thin disc, while the thick-disc slope is null 
within the error. We note that the highest [Mg/Fe] bins are 
more scattered as a result of lower statistics. In contrast, the 
vertical velocity {right panel) shows a first sequence with 
no trend for the thin disc, while a second sequence, which 
corresponds to the thick disc, is characterised by a negative 
gradient equal to —40± 26 km s“^ dex“^. Because we took 
into account the statistics errors to compute the slopes, this 
negative gradient is well explained by the presence of the 
lowest [Mg/Fe] less populated bin. If we derive the slope 
without taking care of the errors on the mean velocity, we 
find a slope equal to —21 ± 30 km s“^ dex“^. Thus, these 
gradients found for the radial and vertical average velocities 
with respect to [Mg/Fe] are consistent with those obtained 
with respect to [Fe/H]. 

The azimuthal average velocity {middle panel) shows a 
strong anti-correlation with respect to [Mg/Fe] that is in 
fact related to the increase of the velocity dispersion with 
age (see Sect. |^. Moreover, the thin and thick discs show 
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two distinct sequences. The thin disc is characterised by 
a negative gradient of —63 ± 10 km s”'^ dex“^ (Spear¬ 
man’s rank correlation coefficient equal to -0.73). Here 
again, we tried to quantify the impact of the [Fe/H]- 
poor and [Mg/Fe]-rich tail of the thin disc. When we 
removed this population, we found a gradient equal to 
—61 ± 11km s“^ dex“^, not significantly different from 
the full thin disc. The thick-disc anti-correlation is twice 
as steep and reaches —172 ± 35km s“^ dex“^ (Spear¬ 
man’s rank correlation coefficient equal to -0.84). This az¬ 
imuthal velocity gradient ag r ees v ery well with the one 
found by [Recio-Blanco et al. ([2014 ) with GES iDRl data 
(VV,^/V[Mg/Fe] = -181 ± SSKTT-i dex"!). 


4.3. Galactic disc velocity dispersion vs. [Mg/Fe] 

We now first present the relation between the Galactic ve¬ 
locity dispersion and the [Mg/Fe] ratio without considering 
the separation between the thin and thick discs to explore 
a possible continuity between the two components in terms 
of velocity dispersion. The resulting relation a vs [Mg/Fe] 
is shown in Fig. Again, the different [Fe/H] curves are 
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Fig. 6. Average radial, azimuthal, and vertical velocity (< Vr >, < V,f, >, and < Vz >) as a function of the [Mg/Fe] ratio. The 
[Fe/H] curves are colour-coded, and the average values of each one are written in the legend. The symbols, errors bars, and linear 
fit are the same as in Fig. 


colour-coded, while the overall trend obtained considering 
all the [Fe/H] range is represented in black. 

For the radial velocity dispersion an, two regimes are 
observed. First, increases for all [Fe/H] curves. The over¬ 
all trend (in black) reaches from 21 ±3 to 59 ±2 km s”^ be¬ 
tween [Mg/Fe] = —0.13 to -F0.28dex. Second, for the most 
[Mg/Fe]-rich bins, the velocity dispersion seems to decrease. 
For the overall trend, a decrease of 9 km s”'^ is observed 
within 1 — a. For the curve with [Fe/H] = —0.15 dex (in 
red), CTfl clearly decreases by 13km s”'^ within 2 —cr. While 
for the curve with [Fe/H] = —0.53 dex (in blue) the de¬ 
crease is less obvious, ctr shows a strong drop for the most 
[Fe/H]-poor curve, 23km s”^ within 2 — cr. We note that 
all the [Mg/Fe] bins showing a decrease are populated by 
more than 19 stars. 

For the azimuthal velocity dispersion cr^, a single be¬ 
haviour seems to be observed. A smooth increase of the 
velocity dispersion is seen for the overall trend (in black) 
from 17±2 to 33±3km s”'^ throughout the [Mg/Fe] range. 
The decreasing velocity dispersion of the highest [Mg/Fe] 
bins observed for is less evident here. A decrease by 
9km s”'^ is detected for the last [Mg/Fe] sub-sample of 
the curve with [Fe/H] = —0.34dex, within 2 — cr. There are 
104 stars. In addition, a drop of is suggested for the two 
most [Fe/H]-rich curves (red and magenta), but it might 
equally well be interpreted as a plateau. 

Finally, we measure a robust increase of az throughout 
the [Mg/Fe] range. The overall trend (in black) reaches from 
13±2 to 37±3km s”'^. The different curves of [Fe/H] follow 
this main trend without showing any decrease as a function 
of [Mg/Fe]. 

For the vertical velocity component, several quantities 
are plotted in Fig. In panel a), the measured velocity 
dispersion is shown as a function of the [Mg/Fe] ratio. For 
each bin, we computed a simple standard deviation of the 
sample, as presented in panel b). Both measured and intrin¬ 
sic dispersions are compared in panel c), and as expected, 
fmeas Is lower. The behaviour of the error on (Tint is also 
shown in panel d). The errors are much smaller than the 
dispersion. For a given curve of [Fe/H], the I/—shape of the 
curve is due to the smaller statistics at each edge. 

To check the reliability of all these velocity dispersion 
trends, we tested the effect of a slight shift by ±0.05 dex on 


the [Fe/H] curves. For both increasing and decreasing 
behaviours are recovered as a function of the [Mg/Fe] ratio. 
For the azimuthal component, we still measure an increase 
of (70 , and the dispersion seems to reach a plateau for the 
highest bins in [Mg/Fe]. Finally, the smooth rise of az is 
recovered on average. 

In addition, we checked the properties of the stars com¬ 
posing the bins that show a decrease ([Mg/Fe] > 0.2dex). 
On one hand, the median errors in Tea and therefore in 
distances are higher than the typical thin-disc errors, but 
still lower than 55 K and 0.13kpc, respectively. On the 
other hand, the median internal errors in log(g), [Fe/H], 
H, and Galactocentric velocities are totally comparable 
with the rest of the sample, and no peculiar lines-of-sight 
were chosen to observe these stars. We also checked the 
median internal error for each [Mg/Fe] bin, which is 
around 0.04 dex, while the most Fe-rich/Mg-poor bin 
reaches 0.07 dex. Moreover, the standard error of the mean 
[Mg/Fe] value for each bin is small, typically 0.006 dex, 
reaching 0.03 dex for the most Fe-poor curve. As a conse¬ 
quence, the velocity dispersion trends are not sensitive to 
the errors on [Mg/Fe]. Finally, as developed in Sect.l^ the 
velocity measurements for low-statistics samples do not 
reflect any bias leading to an artificial under-estimation 
because a minimal number of 15 stars per bin was fixed. 

Considering the thin- to thick-disc separation deter¬ 
mined in Sect.[^ we computed the new velocity dispersions 
for each velocity component of each disc to see whether 
the chemical separation can highlight a possible dynamical 
distinction. The results are shown in Fig. [^ Compared to 
the previous plots without separation, we clearly see that 
regardless of the velocity component, the results are very 
similar. We note that for a given [Fe/H] curve, a [Mg/Fe] bin 
located at the thin- to thick-disc transition is divided into 
two smaller bins, with two values of the velocity dispersion 
on average equal to the entire bin (albeit with larger errors). 
As a result, we can determine that the thin- and thick-disc 
populations show very different behaviours, at least for ct/j. 
The thin-disc sequence follows the main increasing trend of 
CTfl, while the highest values of a and the decreasing bins of 
[Mg/Fe] belong to the thick-disc stars. For and cr^ the 
thick discs can be seen as an extension of the thin disc. The 
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Fig. 7. Radial, azimuthal, and vertical velocity dispersions (crn, , and az) as a function of [Mg/Fe] ratio. The colour code and 
legend are the same as Fig. The errors bars come from the standard deviation of 1000 bootstrapping realisations (as explained 
in Sect. [^. 






Fig. 8. a) Measured vertical velocity dispersion v.s. [Mg/Fe]. The error bars correspond to the standard error of the standard 
deviation, b) Intrinsic vertical velocity dispersion v.s. [Mg/Fe], with error bars from the bootstrap realisations, c) Measured v.s. 
intrinsic vertical velocity dispersion, d) Errors of the intrinsic vertical velocity dispersion v.s. [Mg/Fe]. 





Fig. 9. Same as Fig. but adopting the thin- to thick-disc separation as determined in Sect. (open squares and hlled circles, 
respectively). 


velocity increases smoothly from the thin- to the thick-disc 
sequences without a major visible break. 

As the volume covered by our data extends to distances 
quite far from the Galactic plane (jZj > 1.5kpc), we tested 
a possible contamination of our sample by halo stars that 
could have artificially increased the mean velocity disper¬ 
sion of some metal-poor bins. To do so, we computed the 
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velocity dispersion for three ranges of Z. First, we selected 
the stars with jZj < 0.5kpc (80% of the sample). The typ¬ 
ical increase of a is observed for all the velocity compo¬ 
nents, while the decreasing behaviour is not seen because 
of insufficient statistics in the higher [Mg/Fe] bins. Then, 
we selected the stars with jZj < 1.0kpc (94% of the sam¬ 
ple), recovering most of the same trends as shown in Fig. 
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thanks to the increasing number of high-[Mg/Fe] stars. Fi¬ 
nally, computing the velocity dispersion for the stars with 
jZj < 1.5kpc (more than 99% of the sample), we find the 
trends obtained Fig.[^ as expected. Therefore, the contam¬ 
ination from halo stars does not seem to be responsible for 
the observed patterns that seem to be present close to the 
plane, although with lower thick-disc star statistics. In ad¬ 
dition, to bias the mean velocity dispersion estimations, the 
contamination would need to be in-homogeneous in [Mg/Fe] 
for a given [Fe/H] curve (more halo stars contaminating 
the intermediate [Mg/Fe] bins than the high [Mg/Fe] bins 
presenting lower dispersion values), which seems unlikely. 
Finally, we emphasize that the selection function is not op¬ 
timised to target very metal-poor stars in large amounts, as 
is visible from the [Fe/H] coverage of Fig. Hand Fig. We 
therefore conclude that the contaminationby halo stars is 
not significant in our study. 

In this section, we thus derived robust trends between 
the Galactocentric velocities with respect to [Fe/H] and 
[Mg/Fe] a nd fo un d a g ood consistency with previous stud¬ 
ies (Sects. [tT] and 4.2). Moreover, we established a robust 
relation between [Mg /Fe] and the Galactocentric velocity 
dispersion (Sect. 4.3). This relation is derived here for the 


first time for a very broad spatial coverage, probing both 
above and below the Galactic plane and including the So¬ 
lar neighbourhood. We also showed that the stars with de¬ 
creasing velocity dispersion belong entirely to the thick-disc 
sequence. 

5. Effect of the spatial coverage 

In this section, we first compare the [Mg/ Fe] — a relation 
establis hed above with the recent results by jMinchev et al(] 
(2014b I for the RAVE survey. Then we study in a second 
part the impact of the R — Z coverage on our [Mg/Fe] — a 
relation. 

5.1. Comparison with the RAVE survey 



R ~ 18000. Furthermore, RAVE established the relation 
between the velocity dispersion and the [Mg/Fe] ratio for gi¬ 
ants (log (g) < 3.5) with a ty pical uncertainty of 0.15 dex in 
[Mg/Fe] ( Boeche et al [201 1 1 , larger than the ones adoptee 


in the present work based on higher resolution spectra. 
Minchev and collaborators focused on the Solar suburb, se¬ 
lecting stars in the range \Z\ < 0.6kpc and 7 < R < 9kpc. 
We emphasize that our study covers a larger Galactic vol- 



crease of as a function of [Mg/Fe], regardless of the 

velocity component for [Mg/Fe] < -1-0.1 dex, followed by a 
steeper increase of up to [Mg/Fe] < -1-0.3 dex. We 

also detect a decreasing behaviour for the most [Fe/H]- 
poor and [Mg/Fe]-rich stars for the radial component. How¬ 
ever, we are able to measure a clear decrease of or for the 
curve with [Fe/H] = —0.15 dex because of the hi gher pre¬ 
cision i n [Mg/ Fe]. On the other hand, contrary to Minchev 
et al. (2014bI, we do not detect any decrease of with 


[Mg/Fe], finding a smooth continuous increasing trend. We 
note that the lack of GES stars with [Mg/Fe] > -1-0.4 dex 
m ay be responsible for no t resolving the decrease observed 
by jMinchev et ah (2014bI in this range of [Mg/Fe]. 

To complete our comparison, we now derive a new 
(T — [Mg/Fe] relation, but selecting stars with the same 


Minchev et al. ( 2014b[ | (see the 


R — Z coverage as in 
boxes in Fig. ffl. We kept exactly the same mapping of the 
[Mg/Fe] — [F^H] plane, as well as the same conditions on 
the Galactocentric velocity errors. For each velocity compo¬ 
nent, the number of stars was then reduced, corresponding 
to 48, 54, and 44 % of the three samples presented at the 
end of Sect. [^ We show our results in Fig. |T0| First, we 
see that the stellar populations with low [Fe/H] and high 
[Mg/Fe] tend to disappear. There is no statistics to derive 
a robust velocity dispersion value, the bins sometimes con¬ 
tain no more than three to four stars. Second, the overall 
trend of a increases smoothly with [Mg/Fe] for the three 
velocity components, as expected. Moreover, cjr seems to 
decrease for the highest [Mg/Fe] bin of the grey curve with 
[Fe/H] = —0.35 dex. Despite the lower number of bins avail¬ 
able to analyse with this sub-sample of the data, we see 
that the properties of the velocity dispersion observed with 
RAVE are generally compatible with the iDR2 GES data 
over the same Galactic volume. The GES data allow us to 
extend these behaviours to a much larger volume. 


5.2. Probing the R — Z space 

The broad spatial coverage of the GES iDR2 data allows 
analysing the possible dependencies of the velocity disper¬ 
sion on the position of the stars in the Galaxy (i? and Z). 

For this purpose, we first computed the velocity disper¬ 
sion of each Galactic velocity component as a function of the 
Galactic radius R, taking into account the chemical thin- to 
thick-disc separation (see Sect.[^. In addition, following the 
results of Sect. |4.3 [ 
samples for [Fe/H 


we divided the thick disc into two sub- 
< —0.45 dex and [Fe/H] > —0.45 dex. 
We show our results in Fig. [m First, the behaviour of ctr 
with R for the thin and thick discs shows no clear trend; 
it is nearly constant around 35 and 45 km s“^, respectively. 
For the thick disc, we observe that the velocity dispersion of 
the [Fe/H]-poor population is higher than the [Fe/H]-rich 
population. Second, shows a clear decrease with R for 
both discs. Within 2 — cr, the thin disc is characterised by a 
decline of ~ 15 km s“^ over the whole range of i?, while the 
thick disc seems to decrease by ~ 18 km s“^ in the same 
range. For cr^, no particular trend seems to emerge for ei¬ 
ther the thin or the thick disc. We conclude that the most 
significant trends belong to the thin and the thick disc in 
the azimuthal component with a decrease of with in¬ 
creasing Galactocentric radius. We also note that the trend 
in tT 0 and az may be partially hidden because the sam¬ 
ples in V^ and V^ do not extend as far in R as the one in 
Vr extends more towards the centre and lacks the anti¬ 
centre, while V^ is restricted to 7 —9kpc). However, a clear 
decrease is seen for the thick-disc azimuthal velocity with 
increasing radius (see Fig. 1^). We also find a rough esti¬ 
mate of 6.5 kpc for the radial scale length of the azimuthal 
velocity dis persion in the whol e thick disc. This estimate 
agrees with Lewis & Freemanj (|l989 1 = 6.7kpc). 

We note that ctz is expected to deceas e with R (see for 
example the study of Bovy et al. ( 2012b[ | with the SEGUE 
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[Mg/Fe] [Mg/Fe] 


Fig. 10. Same as Fig. [t] but the spatial coverage corresponds to the one adopted by|Minchev et al.| (|2014b|) for the Solar suburb. 
We note that in this case, the number of stars in our samples is a factor of 4-5 smaller than the one used by|Idinchev et al.| (|2014b|), 
which results in not resolving the decline in az- 


data); we do not observe this behaviour here. We checked 
the case where we considered the disc as a whole (no thin- 
to thick-disc separation) and measured a constant trend of 
az with R. From the trends Fig. El we measured the ratio 
apjaz for both the thin and the thick disc at i? = 8kpc, 
finding apjoz ~ 0.7 for both discs. Our sample being com¬ 
posed of cool d warfs and giants, th is value agrees well with 
o'rIo'z 0.6 of Binney et al. (2014) with the RAVE survey, 
considering no thin- to thick-disc separation. 

We then determined the behaviour of the velocity dis¬ 
persion with the distance to the Galactic plane. To avoid 
possible observational biases, for example the radial trend 
of a presented in the previous paragraph, we selected the 
stars in the range 7 < R < 9kpc. The results are shown in 
Fig. The thick-disc stars are divided into the same two 
sub-samples as in Fig. |ll| First, from a global point of view, 
we did not detect any significant trend with jZj for the thin 
disc (black curve). While ap seems to decrease with Z, a^ 
and az tend to show a higher velocity dispersion at larger 
Z. Second, we found that the velocity dispersion of the thick 
disc is higher than that of the thin disc, regardless of the dis¬ 
tance to the Galactic plane. The complete thick-disc sample 
(orange curve) does not show any trend for ap. In addition, 
az is quite constant (isothermal) for the thin disc, which 
could be the signature of stars with nearly circular orbits. 
A decrease of 10 km s“^ within 2 — cr is observed for a^ in 
the range of Z covered. On the other hand, the trends fol¬ 
lowed by the two thick-disc sub-samples are more scattered 
in the radial component. We note the high-velocity disper¬ 
sion of the thick-disc bins at \Z\ < 0.3 kpc. For a^f, and az, 
the two thick-disc populations are very well separated with 
no particular trends. We thus conclude that the most sig¬ 
nificant trend is the possible decrease of a^ with Z for the 
thick disc, while the thin disc does not seem to follow any 
trend, apart from the higher values of the velocity disper¬ 
sion at large Z. We checked the corresponding values of the 
thick-disc azimuthal velocity for the studied bins and found 
a weak gradient VV^/V^ = —10 ± 7 km s~^ kpc ~^ within 
the 1 — cr error (see Fig. [T^ . 

VV^\ 


Kordopatis et al. (20111 found 
kpc“^, 


a negative gradient VV^VZ = — 19 ± 8 kin s 
but the range in Z was four times broader. A significant 
decrease of Vai with 

eFlRldMIl 


restri cted selection of stars, only for \Z\ > 1.5kpc. Binney 
et al. (20141 observed a strong increase of ap and az with Z 
in the RAVE data, and we do not observe such a strong in¬ 
crease in our data, probably because we separated the disc 
into two thin- and thick-disc components. As a test, we 
checked the behaviour of ap^^^z with Z without separating 
it into thin- and thick-disc parts. On one hand, we observe 
a constant behaviour of ap with Z around 34 km s“^. On 
the other hand, a^ and az increase over the whole range 
of Z of 6 and 11km s“^, respectively. The increasing trend 
of az with Z ther efore agrees with the observations of |Bin-| 
ney et al. (2014) with the RAVE data and also with the 
chemo-dynamical model of Minchev et al. (2013). 


6. Implications on the Galactic disc formation 

In this section, our results are discussed in the context of 
the different Galactic disc formation scenarios. First of all, 
our analysis can be summarized in the following points: 

- We chemically identified two distinct populations in the 
[Mg/Fe] us. [Fe/H] plane and defined them as the thin 
and thick disc (see Sect.[^. 

- As detailed in Sects. [TT] and [T^ we confirm that there 
is no correlation between Yp and Yz with respect to 
[Fe/H] or to [Mg/Fe] for either the thin or thick disc. 
However, we found a weak positive gradient of with 
[Fe/H] for the thin disc, while the thick disc shows a 
strong positive trend. Moreover, is strongly anti¬ 
correlated with [Mg/Fe] for both populations. 

- The thin-disc velocity dispersion increases with [Mg /Fe] 
and seems to display a single regime (see Sect. 4.3). 

- The thick-disc velocity dispersion seems to show a sin¬ 
gle decreasing trend with increa sing [Mg/Fe] at least for 
the radial component (see Sect. |4.3| . The decline is ob¬ 
served for three curves in [Fe/H] for cTij (between 10 and 
20km s“^), while it is only suggested for a^. No decline 
is observed in az- 

- For the thick disc, a^ decreases as a function of R by 
18km s“^, while a^ seems to decline (10km s“^) with 


Z (see Sect. 5.2). 


Z was also observed by Recio-Blanco 
with GES iDRl data for a different and less 


The thin- and thick-di sc dichotomy is still a matter of 
debate (Bovy et al.||2012'a). On one hand, the similarity in 
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Fig. 11. Radial, azimuthal, and vertical velocity dispersion ((t_r, , and az) as a function of the Galactocentric radius R. The 

thin and the thick discs are represented by open squares and filled circles, respectively. The full orange line corresponds to the 
whole thick disc, while the red and blue dashed lines represent its decomposition into two sub-samples (Fe-rich and Fe-poor). The 
error bars correspond to the standard deviation of 1000 bootstrapping realisations. The average [Fe/H] and [Mg/Fe] values of each 
population for the whole R range are marked. 
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Fig. 12. Radial, azimuthal and vertical velocity dispersion {cru, acj, and az) as a function of the height to the Galactic plane. The 
symbols, colours, and error bars are the same as in Fig. EH 


the physical properties between the [Fe/H]-rich thick disc 
and the thin disc around [Fe/H] ~ — 0.15dex suggests a 
common evolution scenario at a certain epoch in the history 
of the two components. Th is b ehaviour was also observ ed by 
Recio-Blanco et al. (20141 and Adibekyan et al. (20111 from 
the beha viour of the la/Fe] abu ndances. Moreover, a recent 
study by Nidever et al. (2014) with the APOGEE survey 
separated the Galactic disc into an a-rich and an a-poor 
population and highlighted an overlap between these pop¬ 
ulations at [Fe/H] 4-0.2 dex. In our study, this common 
behaviour is reinforced by the similar azimuthal and ver¬ 
tical velocity dispersions in the metal-rich regime between 
the two sequences. On the other hand, the chemical distinc¬ 
tion of the two populations in the [Mg/Fe] vs. [Fe/H] plane 
for [Fe/H] < —0.20dex suggests two different evolution¬ 
ary paths in the formation of the Galactic disc. Moreover, 
we clearly detected two regimes in our chemo-kinematical 
analysis, both in the study of the Galactocentric velocity 
trends with [Fe/H] and [Mg/Fe]. The anti-correlation 
of the [Mg/Fe]-low part of the thin disc with [Fe/H] may 
be considered as a consequence of an outward migration of 
these stars due to epicyclic motions with their birth radii 


(Loebman et al.|2011 1. The [Mg/Fe]-rich tail azimuthal ve- 
locity of the thin disc shows no trend with V^, suggesting 
a different evolutionary path than for the rest of the thin 
disc. In addition, both the thin and the thick discs show 
an anti-correlation of with [Mg/Fe], with a break at 
[Mg/Fe] = 0.2 dex. As the velocity dispersion is expected 
to increase with age, older stars ar e expected to show low er 
(see for example Fig. 16 from Haywood et al. (20131). 
Here, the [Mg/Fe] plays the role of the age, and both gradi¬ 
ent are due to the driven effect of the age-velocity dispersion 
relation. The knee is due to the break in the age—[Mg/Fe] 
relation (as shown by Haywood et al. (20131). 


In terms of velocity dispersion, while the thin disc seems 
to have slightly increasing trend with [Mg/Fe], the thick 
disc shows a higher dispersion with a clear decrea se for the 
higher [Mg/Fe] values at different [Fe/H] values. Minchev 


et al. (2014bJ) first detected this unexpected decrease of the 
velocity dispersion with RAVE data. They proposed that 
several mergers of decreasing intensity heated the outer disc 
(which would correspond to the highest values of velocity 
dispersion). These mergers led to subsequent radial migra¬ 
tion of stars from the inner disc, resulting from the induced 
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Fig. 13. Top: Mean azimuthal velocity < > as a function of 

the Galactic radius. Bottom: < > as a function of the height 

to the Galactic plane. The symbols and colours are the same as 
in Fig. El The error bars correspond to the standard error on 
the mean velocity. 


spiral structure. Stars from the inner disc migrating out- 
w ards would show the d etected cooler kinematics (see Fig. 2 
. Minchev et al.|2014a I. We note that this scenario is valid 


only if the migration did not have accompanying kinematic 
heating and that stars on more circular orbits preferentially 
migrate. In agreement with this scenario, the stars from 
our study that show this decline do not seem to be pref¬ 
erentially located in the inner (i? < 8 kpc) or in the outer 
{R > 8 kpc) parts of the Galactic disc, but are homoge¬ 
neously distributed as the rest of the sample in terms of line 
of sight. Moreover, the low ap observed for the sub-sample 
with ([Fe/H]) = —0.15dex and [Mg/Fe] = 4-0.27dex is sig¬ 
nificant in our data, correspon ding to the metal-rich thick 
disc. According to Fig. 4 from Minchev et al. (2014bI, the 
age of these stars with [Mg/Fe] « 0.27 dex would be around 
6 Gyr with a birth radius smaller than their actual positions 
(according to their simulation). The larger metallicity cov¬ 
erage of the decrease presented in our data was neverthe- 
less alre ady predicted by the simulations of [Minchev et ah] 
(2014b|| (see their Fig. 3). We point out that the decrease 
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of az at [Mg/Fe] > 4-0.4 dex in their study is not observed 
in our GES data. 

We note that the method of mapping the 2D abundance 
plane developed in both our study and in that of Minchev et 
al. is sim ilar to the approach of "mono-abundance popula¬ 
tions" by Bovy et al. (2012b). They found that az increases 
for SEGUE dwarf stars with high [Mg/Fe] and low [Fe/H], 
but no decrease fo r these stars is observed. Moreover, [Liu fcj 
in de Ven (2012) measured an increase of az with [Fe/H] 
id [Mg/h'e], but the relation is less clear for [Mg/Fe]—rich 


stars. 

On the other hand, the fact that the older stars 
([Fe/H] < —0.15dex and [Mg/Fe] > 4-0.20dex) show a 
low velocity dispersion, comparable to that of the youngest 
thin-disc stars (Fig. left panel), indicates that these ob¬ 
jects were born in a cool environment, that is to say, were 
shielded from strong merger perturbations in the inner disc, 
where the density is high, and were not heated over the 
subsequent 10 Gyr. Moreover, we observe the velocity dis¬ 
persion decline for different [Fe/H] bins, reflecting different 
ages and stellar generations. This fact traces the chemical 
enrichment in iron of the disc. 

We can also compa re our study with the recent work of 
Haywood et al. (2013), who investigated the age structure 
of the stellar populations in the close Solar neighbourhood. 
Their scenario is based on a thick-disc formation through 
a turbulent phase marked by an intense star formation. No 
decrease of the velocity dispersion was observed by these 
authors for the more [Fe/H]-poor stars (the oldest) of the 
thick disc, although we stress that their statistics is much 
smaller than that of our study. One of their main conclu¬ 
sions was_that_th£radia^Mr^ion, in the sense of churning 
(e.g. Schonrich & Binney 2009) did not play a fundamen- 
tal role in the thick-disc formation. Haywood et al. (2013) 
also argued that the thick disc rernained conhned to the 
inner disc. The peak velocity dispersion that we observe for 
[Mg/Fe] 0.3 dex is compatible with an initial turbulent 
disc with higher random velocities. Nevertheless, this pic¬ 
ture has to somehow leave room for the existence of a small 
proportion of thick disc stars with cool kinematics (those 
in the high-[Mg/Fe] bins with the observed dispersion de¬ 
crease) . 

On the other hand, the cool kinematics of the Mg-rich 
tail of the thick-disc stars does n ot support the idea of ac¬ 
creted satellites, as proposed by Abadi et al. (2003). The 
velocity dispersion of the older stars in our study is too 
low compared to the expected behaviour of the thick-disc 
stars in their simulation (> 80km s“^). The presence of a 
cool thick-disc tail could no t be interpreted by the forma¬ 
tion mechanism proposed by Bournaud et al. (2009) either, 
who evoked a clumpy star formation. 'The proposed pertur¬ 
bations in the disc lead to a scatter of both the stars and 
the gas, increasing the velocity dispersions without room 
for the existence of stars with cool kinematics in the higher 
[M g/Fe] bins at each [Fe/H] ratio. The scenario suggested 
by Brook et al. (2004|, who argued that a gas-rich merger 
was responsible of the thick-disc formation, cannot explain 
the presence of such dynamically cool stars either. In their 
numerical simulation, the thick-disc phase is characterised 
by a strong increase of the stellar velocity dispersion with 
age, again in contradiction with the trends observed here. 

Finally, we can discu ss our results in th e cont ext of 
the scenario proposed by Villalobos & Helmi (2008), who 
argued that the thick disc was built by a pre-existing 
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disc that was heated by a minor merger. Their model 
predicts a remnant-disc with cool kinematics, composed 
of old stars, with a typical location below jZj = Ikpc. 
These old stars would roughly correspond in our study to 
the bins with [Fe/H] < —O.Sdex and [Mg/Fe] > 0.3 dex. 
For a given [Fe/H] bin, these stars are characterised by 
cooler kinematics, and 80% of the stars are located below 
jZj — Ikpc. These stars might be good candidates to rep¬ 
resent this remnant-primary old disc, because it is not af¬ 
fected by the merger event. However, the fact that the stars 
with cooler kinematics show a particular chemical pattern 
(higher [Mg/Fe] values at each [Fe/H] bin) than the stars 
with hot kinematics presents a difficulty for this scenario. 
On the other hand, the possibility of a chemical segrega¬ 
tion between the heated and the non-heated stars could be 
considered. In particular, the old confined disc might have 
been more efficient in forming stars and recycling metals 
than the heated disc. 


7. Conclusion 


From the internal data release 2 of the Gaia-ESO Survey, 
we selected a sample of 7 800 FGK stars to study the chemo- 
kinematical structure of the Galactic disc. We derived the 
kinematics of these stars based on the fundamental pa¬ 
rameters and proper motions. Based on the [Fe/H] and 
[Mg/Fe] abundances, we found a chemical gap between a 
first population with ([Fe/H]) = —0.15 dex and ([Mg/Fe]) = 
-|-0.04dex and a second group with ([Fe/H]) = —0.45dex 
and ([Mg/Fe]) = -1-0.29dex that we associated with the 
thin- and thick-disc sequences, respectively. 

We derived the trends between the three-dimensional 
Galactic velocities and the [Mg/Fe] and [Fe/H] ratios. We 
found a weak correlation between and [Fe/H] for the 
thin disc due to the presence of a [Fe/H]-poor tail, while 
a strong correlation was established for the thick disc, in 
agreement with previous studies. 

Located within 6 < H < lOkpc and jZj < 2kpc, the 
GES iDR2 data allowed us to derive a relation between the 
Galactic velocity dispersion and the [Mg/Fe] ratio, covering 
therefore a large Galactic disc volume. In terms of velocity 
dispersion, on one hand, the thin disc shows a smooth in¬ 
crease of cr with increasing [Mg/Fe]. On the other hand, the 
thick disc can be seen as a continuity of the thin disc and 
reflects a decrease of its radial velocity dispersion for the 
stars with the highest values of [Mg/Fe]. This is the second 
time th at this behaviour is obs erved; it was first discov¬ 
ered by Minchev et al. (2014b| with RAVE and SEGUE 
data. The low dispersion observed for the most [Mg/Fe]- 
rich bins suggests that these stars were formed in a rela¬ 
tively quiet environment. Mor eover, our observations could 
be explained by the scenario of Minchev et al. (2014b): a se¬ 
ries of merger events perturbing and heating the outer disc, 
inducing spiral arms that propel stars outwards in the disc. 
On the other hand, other proposed scenarios of thick-disc 
formation, including external or internal mechanisms, need 
to be refined to take into account this new observational 
constraint: the existence of a sub-sample of old disc stars 
with cool kinematics that show a particular and distinct 
chemical pattern, but are embedded in a generally turbu¬ 
lent context. 

Finally, we conclude that combining individual chemical 
abundances and kinematics is a powerful tool for studying 
the stellar populations of the Milky Way, provided that 


both the chemical and kinematical estimates are derived 
with high enough precision. On the other hand, we worked 
in the context of the [Mg/Fe] used as an age proxy, and we 
note that we cannot exclude that if the [Mg/Fe]-age corre¬ 
lation fails, the stars presenting cool kinematics and high- 
[Mg/Fe] values could be formed at recent epochs, which 
would completely change the picture of the interpretation 
scenarios. Further investigations are necessary to probe an 
always larger volume with an increasing number of stars, for 
example with the next data releases of the Gaia-ESO Sur¬ 
vey and the spatial mission Gaia from the European Space 
Agency. The possibility of deriving precise stellar ages to di¬ 
rectly study the age-velocity dispersion will be a huge qual¬ 
itative improvement in our comprehension of the Galactic 
disc evolution. 
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